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ABSTRACT 

We discuss fragmentation processes which induce star formation in dense walls of ex- 
panding shells. The influence of the energy input, the ISM scale-height and speed of 
sound in the ambient medium is tested. We formulate the condition for the gravi- 
tational fragmentation of expanding shells: if the total surface density of the disc is 
higher than a certain critical value, shells are unstable. The value of the critical density 
depends on the energy of the shell and the sound speed in the ISM. 
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1 INTRODUCTION 

HI shells (or supershells) are structures identified in the dis- 
tribution of the neutral hydrogen (HI) in the Milky Way 
and in many nearby galaxies. The typical shell consists of a 
rarefied medium surrounded by a dense thin wall, in some 
cases expanding supersonically to the ambient interstellar 
medium (ISM). Dimensions of shells lie in the range from a 
few 10 pc to a few kpc. Energies involved in their creation 
are of the order of 10^'^ — lif^erg. Shells have miscellaneous 
shapes, but many of them are nearly spherical, elliptical or 
cylindrical. Observations of shells are reviewed in Brinks & 
Walter (1998). 

The energy needed for the creation of a shell originates 
either in the combined effect of massive stars in an OB asso- 
ciation (radiation, stellar winds and supernova explosions), 
in the very energetic event connected to the gamma-ray 
burst (a hypernova or merging of compact companions in 
a binary system), or in the infall of a high velocity cloud 
(HVC) to the galactic HI disc. In this paper we consider 
the continuous energy input from massive stars of an OB 
association. 

The density of the swept-up gas in walls of shells is 
higher than the average density of the unperturbed ambient 
medium, which increases the probability of star formation. 
Star formation in dense walls of HI shells is observed in some 
galaxies, e.g. in LMC (Kim et al., 1999), SMC (Stanimirovic 
et al., 1999), IC 2574 (Walter & Brinks, 1999), Sex A (Van 
Dyk et al. 1998), etc. Its significance for the evolution of 
galaxies was studied by Palous et al. (1994) and others. 

There are different mechanisms to trigger star forma- 
tion: (1) a compression of pre-existing clouds, (2) an accu- 
mulation of gas into a shell, (3) cloud collisions (Elmegreen, 
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1998), (4) shell collisions (Chernin et al., 1995). Here we fo- 
cus on the mechanism (2) of the gravitational instability of 
mass accumulated in a thin expanding shell. The growth of 
perturbations was analyzed in the linear approximation by 
Elmegreen (1994), nonlinear terms have been included by 
Wiinsch & Palous (2001). 

Any density perturbation on the shell surface is 
stretched by the expansion, while its self-gravity supports 
its growth. The instantaneous maximum growth rate is 
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where R is the radius of the shell, V is its expansion speed, 
Esh is its column density and Csh is the speed of sound within 
the wall. The perturbation grows only if a; > 0. 

In galaxies with thin discs shells differ from a sphere. To 
discuss star formation triggered in expanding shells without 
an a priori assumption on their shapes, we use 3-dimensional 
simulations. In a numerical code, the condition (0) is used 
and we quantify when and where the expanding shell starts 
to be unstable. This approach was first used in Ehlerova et 
al. (1997, Paper I), here we extend parameter ranges and 
generalize results. We also discuss triggered star formation 
observed in nearby galaxies and its connection to expanding 
shells. 



2 THE THIN SHELL APPROXIMATION 

The energy input from an OB association creates a blast- 
wave, which propagates into the ambient medium (Ostriker 
& McKee 1988; Bisnovatyi-Kogan & Silich 1995). During 
the majority of the evolution its radius is much larger than 
its thickness and the thin shell approximation can be used. 
In this approximation the blastwave is considered to be 
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an expanding infinitesimally thin layer surrounding the hot 
medium inside. 



2.1 The analytical solution 

The analytical solution of the expansion in the thin shell 
approximation was derived by Sedov (1959). In a static, ho- 
mogeneous medium without the external or internal gravita- 
tional field, the blastwave is always spherically symmetric, it 
sweeps the ambient matter and decelerates. Neglecting the 
external pressure, and assuming the continuous energy in- 
put L, the self-similar solution for R, V and Esh is (Castor 
et al., 1975): 
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where n and /x are the volume particle density of the ambient 
medium and the relative weight of one particle (see also 
Paper I or Ehlerova, 2000). 

The formula (|^) for the instantaneous maximum growth 
rate of perturbations shows that at early stages of the evo- 
lution the shell is stable, as the fast expansion stretches all 
perturbations which might appear. Only later, when V is 
small, 7? large and high, the self-gravity starts to play a 
role. At the time tb, when the growth rate uj becomes posi- 
tive for the first time, the shell starts to be unstable and the 
fragmentation process begins. 

Inserting the analytical solution (^ - for R, V and 
Esh to the formula we get relations for the time tb, the 
radius R(tb) and the expansion velocity V{tb): 
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The radius R{tb) is a lower limit to the distance on which 
the fragmentation (and triggered star formation) may hap- 
pen; the expansion velocity V{tb) is the expansion veloc- 
ity at which the fragmentation process begins. Since in the 
subsequent evolution the shell further decelerates it gives 
the upper limit to the random component of the velocity 
of newly created clouds (or stars). The relation (^) shows 
that V{tb) does not depend on the density of the ambient 
medium n. This may be one reason why new clouds are born 
in different environments with similar velocity dispersions. 
Relations (j^ , (^ and ^ are derived under the assump- 



tion of the supersonic motion, i.e. 



o{t) > Ce 
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becomes a sound wave before it starts to be unstable (i.e. 
before the time tb), it will always remain stable. Assuming 
that the expansion velocity V at tb has to be greater than 
or equal to the sound speed in the ambient medium Cext , we 
can derive from the equation (Q) the critical luminosity of 
the energy source Lcrit, for which the expansion velocity V 
at the time tb equals the sound speed Cext'- 
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If the energy input is greater than Lcrit, the shell starts to 
fragment; if the input is smaller, the shell is always stable. 
Lcrit does not depend on the density of the ambient medium 
but it is a strong function of its speed of sound. 

Relations (^), (^ and (Q), which apply to the static and 
homogeneous ambient medium, can be extended to spheri- 
cally symmetric cases with the density gradient. Theis et al. 

(1998) show, that the density gradient in the surrounding 
medium has to be shallower than isothermal for the self- 
gravity of the accumulated matter to overcome the stability 
supported by the expansion. 

The analytical solution for shells expanding in the expo- 
nentially stratified ISM was derived by Maciejewski and Cox 

(1999) . It assumes that the shape of the shell can be approxi- 
mated by a prolate ellipsoid. This assumption applies only to 
small shells created by energies corresponding to about one 
supernova in the special ISM distribution. Unfortunately, for 
supershells with energies in the range 

1052 _ 10^^ erg, the 

assumption of ellipticity fails. 

The linear perturbation theory does not estimate the 
evolution of the instability after the time tb very well. A 
better estimate of the evolution a.t t > tb is provided by 
the analysis including nonlinear terms (Wiinsch & Palous, 
2001). 



2.2 Numerical simulations 

The thin shell approximation has been applied in numerical 
simulations in one and two dimensions by many authors 
(see Ikeuchi, 1998, for a review on this subject). Models 
have been further extended into three dimensions by Palous 
(1990) and Silich et al. (1996). 

The code, which we use in this paper, is an improved 
and extended successor to the code of Palous (1990), de- 
scribed in Paper I and in Efremov et al. (1999). The thin 
shell is divided into a number of elements; a system of equa- 



© 0000 RAS, MNRAS 000, 000-000 



Triggered star formation in expanding shells 3 




Figure 1. Fragmentation of shells in Gaussian discs, Cext = 5 kTns~^. The total input energy Etot is given in each panel. The normal- 
ization density parameter = 0.265 cm~^ (left panel); 0.056 cm""^ (center); 0.031 cm~^ (right panel). Results of numerical simulations 
are shown by symbols connected by solid lines, they show the minimum density no needed for the onset of the fragmentation. Dashed 
lines correspond to the constant gas surface density Egas = 14.5 X 10^" cm~^ (left panel); 3.1 X lO'^" cm~^ (center); 1.7 X 10^" cm^^ 
(right panel). 



tions of motion, mass and energy for each element is solved 
and the fragmentation condition (jl|) is evaluated. 

Numerical simulations include some effects, which are 
neglected in the analytical solution. The most important 
are: 1) the pressure of the ambient medium and radiative 
cooling; 2) the finite, time-limited energy input from an OB 
association (life-time of the source is r, we take r = 15 Myr); 
and perhaps the most obvious one 3) the stratification of the 
ISM in galaxies. 

Due to the disc-like ISM density distribution and to 
the form of the gravitational potential, shapes of expanding 
shells are not spherical, but elongated in the direction of the 
density gradient; the differential rotation distorts shells in 
the plane parallel to the galactic plane. In the case of higher 
input energies and thin discs, a blow-out may occur and a 
fraction of the energy supplied by hot stars escapes to the 
galactic halo. Blow-out decreases the influence of the energy 
source on the densest parts of the shell in the galactic plane. 

Fragmentation properties vary with the position on the 
shell. Typically, for shells growing in a smooth distribution of 
gas with the ^-gradient, a dense ring is created in the region 
of the maximum density. This is the most unstable part of 
the shell, while large lobes in low-density regions are stable. 
In the following we present results for the most unstable 
part of shells. Energy sources, which are not located in the 
symmetry plane of the disc, produce asymmetrical shells. 
However, as the fragmentation takes place in the densest 



regions, it is not influenced very much by the asymmetry of 
the shell. The blow-out effects (and the subsequent decrease 
of the "effective" energy) are enhanced in such a case. 

The gravitational fleld of the galaxy may change the 
velocity fleld of the ISM and the shape of the shell, and thus 
conditions for the fragmentation. Simulations show, that the 
differential rotation has a smaller influence than the density 
stratiflcation. Usually, two unstable regions appear on tips 
of an ellipse. Galaxies with little shear (e.g. dwarfs) do not 
influence the shells even to that degree (see Ehlerova, 2000). 

The ISM in galaxies is not perfectly smooth but con- 
tains many small-scale density fluctuations (see Silich et al, 
1996, for the partial analysis of their influence), it is turbu- 
lent etc. However, coherent shells and bubbles are observed 
even in a rather turbulent ISM. Our numerical experiments 
in a medium with a hierarchy of density fluctuations show, 
that shells are mostly influenced by the large-scale gradi- 
ents in the ISM (with the possible exception of the very 
perturbed regions, where the coherent shock does not form). 



3 THE CRITICAL SURFACE DENSITY 

To study the influence of the 2-stratification, we simulate 
shells in different types of the ISM distribution with different 
disc thicknesses. We flx the velocity dispersion in the shell 
Csh (see eq. ^ to a constant value: c^h = 1 kms~^ . Results 
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for the Gaussian stratification n[z) — noe"""^ are given 
in Fig. 1^. Symbols denoting different Etot show the lowest 
density no needed for the onset of the instability during the 
shell expansion. For a given disc thickness H, a total energy 
Etot [Etot is the total energy delivered by the OB association 
to the shell: Etot = L x t) and a velocity dispersion Cext, 
shells are unstable only if no exceeds the depicted value. For 
lower values shells remain stable. Dashed lines are lines of 
the constant gas surface density in the disc (Ejas). 

As can be seen in Fig. hi, lines of the constant disc 
surface density Egaa approximately separate stable and un- 
stable regions. There are two types of deviations between 
the solid line separating stable and unstable regions and 
Egas = const line: 

• The blow-out effect. It is visible in the middle and 
right panels (i.e. in the case of higher energies) for thin discs: 
a higher density no is needed for the instability than pre- 
dicted by Egas = const. The blow-out enables the leakage 
of the energy to the galactic halo, leading to the decrease of 
the effective energy and pressure pushing the densest parts 
of the shell. 

• The small shell in the thick disc. It is visible in the 
left panel (i.e. in the case of a low energy) and thick discs: 
a higher value of no than predicted by Egas = const line 
is needed for the instability. Low energy shells are generally 
small and in thick discs they evolve in an almost homoge- 
neous medium as they never reach dimensions comparable 
to the thickness of the disc. Consequently, the value of the 
gas surface density of the disc is irrelevant. The critical value 
of no has to be higher than predicted by Egas — const cri- 
terion, since a substantial fraction of the ISM in the disc 
remains untouched by the shell. 

A different behavior is observed in complementary sim- 
ulations with r = oo (long-lasting energy supply). In this 
case shells fragment at a radius R{tb) ~ 1.5H. Inserting this 
condition into eq. (6) we get no oc l/_ff^ for given Csh and 
L. The deviation from this behavior in Fig. 1 is connected 
to the fact that t < ti,. 

Other simulations with different disc profiles (Gaussian, 
exponential, multicomponent) show that the critical value of 
the gas surface density T,crit does not depend on the profile. 
Ecrit depends strongly on Cext'- for higher values of c^xt the 
fragmentation starts at higher values of Egas • We derive the 
fit: 

™-»-"(Tife;)"'(iSFT)" (») 

This relation is illustrated in Fig. ti. 



4 COMPARISON TO OBSERVATIONS 

In this paragraph we study one example of the HI hole with 
star formation on the rim. This is the structure number 35 in 
the galaxy IC 2574 (Walter & Brinks, 1999). Taking the mea- 
sured value of Cext ~ 7 km s~^ and the total energy input 
from stellar winds and supernovae of the central star cluster 
inside this shell, estimated by Stewart & Walter (2000) as 
Etot = 4.1 X 10^^ erg, and using the formula (^) we derive the 
critical gas surface density Ecrit = 1-3 x lO^'^ cm~^ . The gas 
surface density in the vicinity of this shell, 1.8 x 10^^ cm~^ 
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Figure 2. The critical gas surface density Scrit in the {logEtot, 
Cext) plane. Various symbols show T,crit from numerical experi- 
ments: * = Ecrit < 10^0 cm- 2; o = 10^0 < Scrit < 4x10^°; A = 
4 X 10^0 < Ecrit < 10^1; □ = 10^1 < Scrit < 2 X 10^1; o = 
2 X 10'^^ < Ecrit- Contour lines are calculated by the formula (]M 
and labelled by the logarithm of E^rit in lO'^" cm~^. 



Name 


Cext 


y^obs 
crit 


Etot 
10^^ erg 


IC 1613 


7.5 


3 


257 


Ho II 


6.8 


6 


92 


DDO 155 


9.5 


3 


640 



Table 1. Values of the minimum input energy to trigger star for- 
mation on the rim of the shell, Etot, as derived from the observed 
values of Cext and E°^^j. 

(as given by Walter & Brinks 1999), is higher than the de- 
rived value Ecrit, which is consistent with our hypothesis. 

We may derive R{ti,) using formula (6) for a homoge- 
neous medium. For this hole we get R{tb) ~ 500 pc, which 
is quite close to the observed dimension. The disc thick- 
ness H of IC 2574 is 350 pc, which is less than R{ti,), and 
the assumption of homogeneity is not fulfilled. In our opin- 
ion, the condition (9) is more suitable for a discussion of 
the fragmentation process since it takes into account the 
z— stratification. 

A comparison of the HI distribution in IC 2574 to the 
Ha emission (which is very often connected to rims of HI 
shells) suggests, that majority of Ha emission lies in the re- 
gions with Egas > 1.7 X lO^'^ cm~^ . If this is the critical gas 
surface density Ecrit, we can estimate the value of the mini- 
mum energy necessary to produce an unstable shell, which is 
Etot = 3.8 X lO^^erg. It means, that in IC 2475 clusters with 
~ 40 massive stars and more can trigger star formation. 

In Table |l| we show the observed average value of the 
velocity dispersion Cext and E"^!^, the azimuthally averaged 
value of the gas surface density, below which star forma- 
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tion rate strongly declines, for three galaxies selected from 
Hunter et al. (1998). Etot in the last column of Table |l| 
derived from the formula (^), gives the value of the mini- 
mum input energy necessary to create an unstable expand- 
ing shell. Corresponding numbers of massive stars in the star 
cluster differ for different galaxies: in IC 2574 and Ho H less 
numerous clusters than in IC 1613 and DDO 155 are needed. 



5 CONCLUSIONS 

We study the gravitational instability of expanding shells 
using the analytical solution and numerical simulations. The 
influence of the total energy Etot, the velocity dispersion of 
the ISM Cext, the thickness and the profile of the gaseous 
disc are analyzed with the following result: the critical value 
of the disc gas surface density J^crit exists; shells expanding 
in discs with the lower gas surface density are stable, shells 
expanding in discs with the higher density are unstable, can 
fragment and form gaseous clouds and later stars. 

Values of Ecrit for reasonable values of Etot and Cext are 
of the order of (10^" - 10^^) cm , which coincides with the 
value of observed threshold surface densities for star forma- 
tion in galaxies (Kennicutt, 1997, 1998; Hunter et al., 1998). 
This agreement may indicate, that the contribution of star 
formation induced by shells to total star formation is impor- 
tant. 

We should distinguish between two cases: 1) the gas 
surface density Egas in the shell is increased compared to 
the average, which decreases the value of Toomre (1964) or 
Safronov (1960) parameter Q for the local gravitational in- 
stability of the galactic disc; 2) the material accumulated 
in the shell starts to be unstable according to the formula 
(0). We label the first case an enhanced mechanism of spon- 
taneous star formation, the second case describes triggered 
star formation. According to our simulations, the main dif- 
ference lies in the dependence on the velocity dispersion in 
the unperturbed medium c^xt'- the critical density Ecrit is in 
the case of spontaneous star formation directly proportional 
to Cext, but in the case of triggered star formation it depends 
on ctxt- A much sharper dependence of Ecrit on Cext is con- 
nected to the mass accumulation to the shell, which stops 
when the shell decelerates to Cext- 

The very steep dependence of T.crit on c^xt indicates 
the importance of the self-regulating feedback for triggered 
star formation mode. Young stars in OB associations release 
the energy and compress the ambient ISM creating shells. 
In dense walls of shells, star formation may be triggered if 
the disc surface density surpasses a critical value Ecrit. Star 
formation is accompanied by the heating of the ISM, i.e. in- 
creasing Cext- This leads to the increase of Ecrit and a subse- 
quent reduction of triggered star formation rate. The energy 
dissipation and cooling of the ISM decreases Cext and Ecrit 
and increases star formation, closing the self-regulating cy- 
cle of the triggered mode of star formation. 

The heating of the ISM and the increase of Cext related 
to star formation is a localized process operating on scales of 
100 pc—1 kpc in size. Consequently, triggered star formation 
is a more local process than the spontaneous mode. The 
less steep dependence of Ecrit on Cext for spontaneous star 
formation means, that this mode may be also effective in 



regions, where Cext is increased and the triggered mode of 
star formation temporarily stopped. 
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